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We consider the consequences of the absence of Birkhoff 's theorem in theories of modified gravity. 
As an example, we calculate the gravitational force on a test particle due to a spherical mass shell 
in the Dvali-Gabadadze-Porrati model (DGP). We show that unlike in General Relativity, the force 
depends on the mass distribution. In particular, the gravitational force within a spherical mass shell 
depends on the geometric structure of the bulk, and is likely non-zero. 
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I. INTRODUCTION 

An important tool in astrophysics and cosmology is 
the ability to reconstruct an object's mass from mea- 
surements of the gravitational force in its vicinity (say, 
by measuring velocities). Newtonian gravity's inverse- 
square law guarantees that the gravitational flux through 
an enclosed surface is conserved, giving us a leverage on 
the mass within that surface; this is Gauss's law. Even 
though General Relativity (GR) is a highly non-linear 
theory, Birkhoff's Theorem (BT) plays the same role and 
allows us to reconstruct the mass of spherically symmet- 
ric configurations. Two particularly important aspects of 
BT are that: (a) outside a spherically symmetric mass 
distribution, the gravitational potential (i.e. metric) de- 
pends on the distribution of the matter density only at 
second order in the potential, due to the role of binding 
energy as a source of gravity; and (b) a shell of mass 
has no effect on the metric in its interior. Without these 
properties we presumably could not calculate almost any 
gravitational fields without knowing details of the dis- 
tribution of matter all over the Universe. Indeed many 
calculations in gravity would become undoable either in 
principle or practice, and many others would become 
enormously more difficult. In general, modifications to 
GR tend to violate BT Q. 

Nevertheless, there has been considerable interest in 
theories that modify GR and break BT. The first of 
two primary phenomenological motivations for interest 
in such theories is the observation that at scales above 
stellar clusters there is too much gravity to be sourced 
by the observed matter. Thus either there is dark matter 
in such systems, or Newton's Law is not valid in this do- 
main. The latter possibility, known in a surprisingly suc- 
cessful, albeit phenomenological form as Modified New- 
tonian Dynamics (MOND), explains a range of such sys- 
tems. The recent discovery of covariant theories 0, Q 
which yield GR, Newton and MOND in appropriate lim- 
its has further fueled interest in modified gravity as an 
alternative to dark matter. 

There is also interest in modified gravity theories as an 
alternative explanation to dark energy for the observed 
accelerated expansion of the universe. Again the exis- 
tence of reasonable covariant extensions to the Einstein- 
Hilbert (EH) action, such as the DGP model of Dvali, 
Gabadadze and Porati (3] and f{R) theories [1], has been 
particularly important. 



Theories that violate BT should also lose the benefits 
afforded by it. It would be expected that the geometry 
locally will depend on the detailed distribution of mass 
in the universe. If this dependence is sufficient to qual- 
itatively affect the behavior of dynamical systems, then 
we could be in a most difficult situation of having alter- 
native gravity theories in which we are fundamentally, or 
at least practically, unable to make firm predictions for 
a wide variety of important quantities. 

We will not perform an exhaustive calculation of non- 
Birkhoffian behavior in all modified gravity theories, al- 
though that may well be a desirable program of research. 
Instead we will begin that endeavor by considering a par- 
ticular example, DGP. We choose DGP both because it 
is a single theory, with only one free parameter of in- 
terest, and because the calculations are readily doable. 
Other modified theories of gravity in general will not re- 
spect BT, and a similar investigation for other theories 
of gravity is in progress Q , with preliminary results that 
are consistent with our expectations that DGP is a rea- 
sonable exemplar. 

To cut precipitously to the bottom line, by solving the 
DGP equations for concentric spherical shells of dust in a 
background with accelerated expansion, we will demon- 
strate the violations of BT that can be expected in mod- 
ified gravity theories. The violation in theories that seek 
to replace dark matter is likely to render the calculability 
of cosmology beyond the homogeneous background sig- 
nificantly more difficult or wholly impossible. Theories 
which replace dark energy may be safer, although even 
there we should still expect significant effects on observ- 
able scales and should not be complacent. 

The paper is organized as follows. In section |TT] we first 
review the background de-Sitter solution in DGP theory, 
and then solve the metric functions on the brane in the 
presence of a weak source on the brane. We show that 
the equations define the metric only up to a function, 
which we name g{r) . We present or modeling for g{r) in 
section IIII[ and then consider and estimate the severity 
of the deviations from BT. We conclude in section HVl 



II. SPACE TIME GEOMETRY 

DGP introduces an infinite flat 5-dimensional space 
(bulk) in which gravitons move freely, while Standard 
Model fields are confined to the 4-dimensional brane 
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world. The 5-dimensional gravity theory retains the EH 
action, while a 4-dimensional EH action is presumably 
induced by radiative corrections by the matter on the 
brane. The action of the theory (other than a possible 
Gibbons-Hawking term) is Q 

+ / d^x^/^C,^ . (1) 



0. Following them, we write the line clement as 



Here z is the extra dimension, A/5 is the fundamental 
five-dimensional Planck mass, Mp is the observed four- 
dimensional Plank mass, g and R are the 5-dimensional 
metric and Ricci scalar, and i?*-^-' and i?*-'*-' are the induced 
quantities on the brane. The gravity mass scale in the 
bulk is taken to be extremely low, M5 ~ 10~^eV. This 
introduces a new physical scale, the crossover scale, at 
which gravity becomes 5-dimensional, rp = Mp/2M^. 
While Einstein gravity approximately holds on the brane 
on distances shorter than tq, one expects modifications 
to gravity at larger distances due to the graviton's free 
movement in the bulk. Actually, bulk effects leak into the 
brane on even shorter distances than ro, due to non- liner 
interactions of the scalar degree of freedom of gravity 
Taking this into account, the length scale at which 
gravity is modified, known as the Vainshtein radius, is 
= (roTs)^/'^, where is Schwarzschild radius of the 
source. Since <SC "Toj this may open an observational 
window to the existence of higher dimensions Q . 

An attractive feature of DGP is that it has a branch of 
cosmological solutions in which cosmic acceleration oc- 
curs without dark energy [lol | as the Hubble function ap- 
proaches ro at late times. By choosing r^^ ^ Hq, DGP 
presents an alternative mechanism for the present accel- 
eration. We follow and solve in DGP for the metric 
of a spherically symmetric, static matter source in a de- 
Sitter background. We then check how the gravitational 
field changes when we keep the mass constant but change 
its distribution. We show that there may be a gravita- 
tional force within a spherical mass shell, depending on 
the form of an undetermined function, g{r). 

We start with the Einstein equations in the bulk: 



(2) 



Here Gab is the bulk Einstein tensor, G^g is the intrinsic 
Einstein tensor on the brane. Tab is energy momentum 
tensor. We have explicitly specified that Tab vanishes 
except at the location zg of the brane in the extra di- 
mension (see below, Eq. (O). Without loss of generality, 
we can take zq = 0. 

We wish to find the metric (at least on the brane) 
induced by a static, spherically symmetric (3-d) weak 
matter source on the brane, in a background cosmology 
which, given observations, we take to be de Sitter-like. 
This problem was first addressed by Lue and Starkman 



N^{r, z)dt^ - A^{r, z)dr^ - B^{r, z) 
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(3) 



We assume a Z2 symmetry across the brane, and that 
there is no spatial curvature in the brane directions, 
B\z=o = r. The resulting bulk Einstein tensor and Ein- 
stein equations are presented in Q (equations (2.7) and 
(2.8) respectively), and we do not reproduce them here. 



A. Vacuum solution 

The metric functions for the background de-Sitter so- 
lution 

7V(r,z) = {1tH\z\){1~ H\^)^^^ 

A{r,z) = {ItH\z\){1- H\^y^^^ (4) 
B(r,z) = {lTH\z\)r, 

with the vacuum energy-momentum tensor given by 

T^ ^ S{z)dmg{pH, -Ph, -Ph, -Ph,0)- (5) 

The energy density and pressure due to the bulk are 



PH = -PH = ^^{roH±l) 
87rrn ^ ' 



B. Weak source on the brane 



(6) 



Adding a gravitational source on the brane, the energy 
momentum tensor remains of the form ([5|), with 



p = PH+Pa{r), P = PH + Pg{r) 



(7) 



Here Pg{r) and Pgir) are the energy density and pressure 
due to the distribution of matter on the brane. It is 
convenient to redefine the metric functions, 



N{r,z) = l + n{r,z) 
A(r, z) = 1 + a{r, z) 
B{r,z) = {l + b{r,z))r . 



(8) 



It can be seen from Eq. ([4]) that to first order in r and 2, 
we get n H'^r^ H\z\, and the solution for empty space 
has n(r, 0) = ztH (dots denote d/dz). Treating the mass 
distribution on the brane as a perturbation around the 
vacuum solution we can write 



n(r,0) = ^H + g{r) , 



(9) 



where g{r) is the correction due to the mass distribu- 
tion. It is here that we deviate from which took 
g{r) = 0. Even thou gh q (r) is not strictly defined by 
the equations, and [13|| argue that g{r) = is essen- 
tial to recover a well-behaved bulk solution. This has 
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been a matter of some dispute, with some authors 
proposing alternative Schwarzschild solutions. Our own 
calculations suggest that the solution of [§| itself requires 
a small non-zero g{r) to recover the desired bulk solution 
for self-consistency. 

We assume that the mass distribution is of negligible 
pressure, Pg{r) w 0, and we define Rg and Gg to be 
the 3- volume integrals over the mass density and g{r) 
respectively, 



i;2 



g{r)df 



(10) 
(11) 



We now want to solve for the metric functions a, b, and n. 
Except for Gzz and Gzr, the bulk equations can be sat- 
isfied by choosing suitable quadratic terms. It has been 
shown that Gzr is identically zero when covariant con- 
servation of the matter source and the brane boundary 
conditions are considered @. The equations we arc left 
with (shown to first order) are the bulk equation Gzz 

n" + ^-2{^y = {±H-g){a + 2b)-{2a + b)b (12) 

(primes are derivatives with respect to r and dots are 
derivatives with respect to z), and the Einstein equations 
on the brane: 



^0 



i + 2b 
2b 

a + b = ro ( 



,2a 
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/2a 2n' , 
ro(^ )-H{ZHr^±2 
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3H{roH±l) 

9 (13) 
-H{3Hro±2) - g. 



Eliminating d and b and integrating dr yields 

nV^ + ra + 
We define a new quantity. 



Ga — Ra 



. 



(14) 



With this solution for /(r) we get simple expressions for 
the metric functions on the brane: 



rt(r, 0) = TiQ + 



dr 



fir) 



a(r,0) 
&(r,0) 



4ro ^ 2rJ . 

i(i?g(r) - Gg{r)) - ^H\^ - rn'{r, 0) 
. 



(17) 



This solution reduces to that of Q in the limit g{r) ~ 0. 
We have generalized by allowing the influence of the 
bulk to be space-dependent, as one expects. While g{r) 
approaches zero far from the source, Gg{r) reaches a con- 
stant non-zero value. The effect of Gg is to modify Rg{r), 
and the sign of Gg{r) will determine whether the gravita- 
tional field will increase (suggesting a gravitational source 
of greater than actual mass) or decrease. The quantity 
relevant for the acceleration of a test particle is n'(r, 0), 
given by 



n'(r,0) ^ 



r 

9{r)r 
4ro 



— (Rg{r)-Ggir)) 
^ 2ro' 



(18) 



III. 



GRAVITATIONAL EFFECTS OF MASS 
SHELLS 



We shall see that the effects of matter on the metric 
are not confined to the region exterior to it. To see this 
we study the gravitational force exerted by a spherical 
mass shell in both its interior and exterior with a reason- 
able model of g{r). We next superpose several concentric 
shells, and construct a mock model of a galaxy, within a 
cluster and super-cluster for example. 



A. modeling 



— (Rg-Gg) 



and use it to replace b and a-|-26 in (|T2|) and p3|) . Substi- 
tuting into multiplying by and integrating with 
respect to r yields a quadratic equation for f(r). Im- 
posing the boundary condition that asymptotically the 
solution approaches the de Sitter background yields 



fir) 



r 



D{r) = (3-2ro(gTi^))V 



^2 r 



+ 16- 



^-r{?>-2roigTH)) 
1 



(16) 



12gr^ro 



Q+^iRg-Gg)+7-\2H^±^) 

ZTq 



Q{r) EE THiroGg - ^gr') - ^ j\{f)g\f)f^df 



Based on dimensional analysis we expect a shell 
with mass Mi and radius to modify fi by roughly 
^2GMi/rf . By choosing this dimensional form we rec- 
ognize g{r) as a correction to the local Hubble scale due 
to the mass shell. We therefore model g{r) as 



5(r) = -y^2GM,/r3x || 



ri/i 



for r < Tj , , 
.4 . - ^ (19) 
tor r > Tj . 



Exterior to the shell, 9(7-) needs to decay strongly enough 
to ensure that Gg{r) is finite in the limit r 00, which 
is satisfied with our choice of r^^. As can be seen later, 
the dominant contribution of g{r) comes from the interior 
region, and the details of the decay external to the shell 
are less important. We have chosen the sign of g{r) to 
augment the attractive force of gravity with an eye to 
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our interest in modifications of gravity that mimic dark 
matter. 

First consider a test particle external to the shell. For 
distances much larger than the shell radius, the only 
relevant contribution comes from the Gg term of (jl7p . 
Asymptotically Gg 2^j2GMirf /r'^, which depends on 
the crossover scale as well as the shell mass and radius. 
Notice that the dependence of Gg on the location of the 
shell is a positive power; hence for fixed mass the contri- 
bution grows for larger shells. This shows explicitly that 
the gravitational force of spherically symmetric systems 
depends on the distribution of matter. 

Now consider a test particle interior to the spheri- 
cal shell. Newtonian gravity predicts that the particle 
feels no acceleration. For DGP, this depends on the bulk 
structure. Taking g(r) non-zero for r < ri modifies the 
gravitational potential and gives a non-zero acceleration 
within the shell. A local deviation of ri{r,0) from the 
DGP-de Sitter value in the presence of a matter concen- 
tration would seem to be expected if, looking at (O, we 
can think of g{r) as a local change in the effective Hubble 
scale. Overdensities are expected to create such modifi- 
cations. This result is extremely disturbing. It suggests, 
for example, that when calculating gravitational forces 
within a galaxy, one must consider the effects of the clus- 
ter and super-cluster in which it is embedded. 



B. Multiple shells 

Notice that g{r) is non-linear in the mass, so one can- 
not superpose solutions. Concentric shells however, each 
contributing a g{r) correction to the Hubble scale, is a 
solvable system. As a mock model of our extended neigh- 
borhood we model our galaxy, cluster and supercluster 
as concentric spherical dust shells, see Fig. [TJ It is self ev- 
ident that galaxies and clusters are not spherically sym- 
metric and concentric shells, and we have not attempted 
any data fitting. Nonetheless, the spherically symmetric 
calculation does give an order of magnitude estimate as 
to how severe are the deviations between two calculations 
within the same theory. As we show, the deviations are 
not negligible. Indeed, another possible way to model the 
galactic neighborhood is with spheres of mass instead of 
shells. This will give yet another possible solution within 
the same theory. As we consider all of these solutions just 
an estimate to the order of magnitude of the problem, we 
chose to work with the simple shells model. 

We first consider DGP as a dark energy alternative, 
and so take tq — 1/ H . For the radii of the shells we take 
lOV, 3 X lO^pc and 3 x lO^pc, and lO^^^f©, IO^^Mq 
and IO^^Mq for the masses. Fig. [2] shows the gravita- 
tional acceleration as a function of the distance. The 
plot compares the acceleration accounting only for the 
galaxy (solid), for the galaxy and cluster (dashed), and 
for the galaxy, cluster and super-cluster (dotted), to the 
Newtonian acceleration (dot-dashed). Clearly, the accu- 
mulating effect of external shells can be non-negligible. 
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FIG. 1: Mock model of a galactic neighborhood as concentric 
spherical shells of mass. 
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FIG. 2: Gravitational acceleration as a function of radius, g 
(solid) includes only the galaxy as a source, (7 + c (dashed) 
adds the cluster, and g-f c-|-sc (dotted) adds the super-cluster. 
The Ne-wtonian acceleration is (dot-dashed). The crossover 
scale is ro = lO^rg = 10^°pc« l/Jfo- 



The implication that mass shells can alter the gravi- 
tational acceleration interior to them seems particularly 
ominous for modified gravity theories that seek to replace 
dark matter as an explanation for the anomalous dynam- 
ics of galaxies. Although DGP is not a candidate MOND 
theory, we wish to understand, in this calculable theory, 
the magnitude of the effects we might expect from BT 
violations. We now let our galaxy, cluster and superclus- 
ter shells have radii of 3 x 10'*pc, 3 x lO^pc and 3 x lO^pc, 
and masses of 10^^ Mq, 10^^ Mq and lO^^M©, account- 
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FIG. 3: As in Fig.0 but with ro = = IC'^pc. 
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ing only for the visible matter. To mimic dark matter, 
To ~ 2.6 X 10''pc is chosen to boost the acceleration due to 
the gravity so that the circular velocity of of a test parti- 
cle in orbit at rg is 200km/s. As expected, smaller values 
of To cause larger modifications to the gravitational force, 
as shown in Figure [3l The cluster and supercluster keep 
the local acceleration approximately constant from the 
radius of the galaxy shell outward. In addition, the in- 
ternal accelerations to the galaxy shell are substantial. 
Clearly "exterior" mass cannot be ignored. 

IV. CONCLUSIONS 

In this work we have used DGP as an example of the 
sickness associated with violation of Birkhoff's theorem. 
We have generalized the solution of 0, and shown that 
(a) even for spherically symmetric systems, the exact dis- 
tribution of the matter affects the gravitational force ex- 
ternal to the source, and (b) one cannot neglect distant 
gravitational sources, because even interior to a spheri- 
cally symmetric mass shell the gravitational force is non- 
zero. These effects pose severe conceptual and calcula- 
tional difficulties. 

The problem for DGP may be due to our ignorance of 
the bulk. Conceivably a better understanding of how the 
4-dimensional metric of the brane is induced and of bulk 
effects may change our grim conclusion and answer satis- 
factorily the issues we have raised. The full understand- 
ing of the gravi-scalar degree of freedom near a brane 



source should hopefully alleviate the ambiguity in the 
determination of g{r)] as argued for example by and 
[l3l |. the gravi-scalar is frozen out, yielding g{r) = 0, ef- 
fectively enforcing a Bikrhoff-like vanishing of the grav- 
itational acceleration inside a spherical shell. However, 
we may find ourselves with a well-determined solution 
which has a non zero g{r) and thus Birkhoff-violating. 

This sickness is not unique to DGP theory: most theo- 
ries of modified gravity are likely to suffer the same prob- 
lem in various degrees. Since DGP is not a candidate 
MOND theory, it is not primarily in their application to 
DGP per se that our results should be viewed, rather 
as an example of the magnitude of the problem that is 
likely to apply to most modified gravity theories. These 
issues are the subject of further investigation. In future 
work we will also consider deviations from both spherical 
symmetry and smoothness and how they affect the grav- 
ity in modifications to GR. A possible conclusion of our 
work may be that we need Birkhoff's theorem to hold for 
a theory to be calculationally tractable. Nature may, of 
course, refuse to cooperate. 
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